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Equilibrium and out-of-equilibrium transitions of an off-lattice protein model have been 
identified and studied, fn particular, the out-of-equilibrium dynamics of the protein under- 
going mechanical unfolding is investigated, and by using a work fluctuation relation, the 
system free energy landscape is evaluated. Three different structural transitions are iden- 
tified along the unfolding pathways. Furthermore, the reconstruction of the the free and 
potential energy profiles in terms of inherent structure formalism allows us to put in di- 
rect correspondence these transitions with the equilibrium thermal transitions relevant for 
protein folding/unfolding. Through the study of the fluctuations of the protein structure 
at different temperatures, we identify the dynamical transitions, related to configurational 
rearrangements of the protein, which are precursors of the thermal transitions. 
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§1. Introduction 

Biopolymers such as proteins and nucleic acids are paradigmatic examples of 
complex systems. Similarly to glasses and to super-cooled liquids, they are character- 
ized by complex free energy landscapes (FELs), which determine their dynamical and 
thermodynamical properties. Manipulation experiments on single biomolecules have 
made it possible to observe unfolding and refolding trajectories of single proteins, 1 ) 
or RNA molecules. 2 )' 3 ^ From a theoretical point of view, the unfolding and refold- 
ing of biomolecules represent typical stochastic processes where out-of-equilibrium 
single trajectories of microscopic systems can be observed, since the typical exper- 
imental time is much smaller than the typical molecular relaxation time. When 
a thermodynamical system is driven far from equilibrium, classical linear response 
theories and other near-equilibrium approximations generally fail. However, in re- 
cent years some general relations for systems driven far from equilibrium by large 
external perturbation have been obtained. These relations, which are known as fluc- 
tuation relations, 4 ) go beyond linear response theory valid only in the vicinity of the 
equilibrium regime. On the one hand, single molecule experiments represent excel- 
lent test-bed for these results obtained in the field of out-of-equilibrium statistical 
mechanics. 3 ) On the other hand, these relations can be used to characterize the ther- 
modynamical properties of biomolecules, overcoming the intrinsic out-of-equilibrium 
nature of the unfolding experiments. In particular, a fluctuation relation introduced 
by Hummer and Szabo 5 ) can be used to estimate the equilibrium free energy land- 
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scape of a system as a function of an internal coordinate: such a relation has been 
used to evaluate the FEL of model 6 ) -8 ) as well as of real proteins. 9 )' 10 -* 

In two previous papers, 7 ) we have shown that the free energy landscape of 
a model protein can also be evaluated by using the inherent structure (IS) ap- 
proach, a method previously used to characterize the structural-arrest temperature 
in glasses 11 ) and super-cooled liquids. 12 ) The investigation of the IS distributions 
allows us to give an estimate of the energetic barriers separating the native state 
from the completely stretched configuration along the out-of-equilibrium unfolding 
trajectories. Moreover, the thermal energies to overcome these barriers are related 
to three temperatures, which are quite similar to the temperatures usually employed 
to characterize the thermodynamical transitions associated to protein folding. 

Thus, the aim of this paper, is twofold. On the one hand we want to reconstruct 
via out-of-equilibrium measurements the free energy landscape of a model protein, 
by using both a fluctuation relation and IS approach. On the other hand, we want 
to characterize the equilibrium dynamical transitions for the same model protein in- 
duced by temperature variation, and to compare these with the information gathered 
via the out-of-equilibrium mechanical manipulations. 

The paper is organized a follows, in Section [2] we briefly describe the model 
protein used in the present work, and the numerical simulations we perform. In 
Section O we introduce the work fluctuation relation and reconstruct the free en- 
ergy landscape by combining it with out-of-equilibrium unfolding simulations. We 
discuss how the FEL can be evaluated via the IS approach in Section U and make 
a comparison with the results obtained in the previous Section In Section we 
discuss the three characteristic temperatures characterizing the thermal unfolding of 
the protein, and relate them to the structural transition as identified in the Sect. 01 
Section [6] is devoted to the analysis of the dynamical transitions observed for this 
model protein, which appear to be precursors of the thermal unfolding transitions. 
We conclude and summarize our results in Sect. [7J 

§2. The protein model 

The model studied in this paper is a modified version of the 3d off-lattice model 
introduced by Honeycutt-Thirumalai 13 ) and successively generalized by Berry et al. 
to include a harmonic interaction between next-neighboring beads instead of rigid 
bonds. 14 ) The model consists of a chain of L point-like monomers mimicking the 
residues of a polypeptidic chain. For the sake of simplicity, only three types of 
residues are considered: hydrophobic (B), polar (P) and neutral (N) ones. 

The intramolecular potential is composed of four terms: a nearest-neighbor har- 
monic potential, V\, intended to maintain the bond distance almost constant, a 
three-body interaction V2, which accounts for the energy associated to bond angles, 
a four-body interaction V3 corresponding to the dihedral angle potential, and a long- 
range Lennard-Jones (LJ) term, V4, acting on all pairs i, j such that |i — j\ > 2, 
namely 



Vi(r i;i+ i) = a(r iji+ i - r ) 2 , 



(2-1) 
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V 2 (9i) = Acos(0i) + Bcos{26i) - V , (2-2) 
V 3 (^, O i} 9 i+1 ) = d[l- S{9 u 9 i+X ) cos(^))] + A[l - S^i, m ) cos(3^))], (2-3) 

tt(n*) = ^(4-3M • (2-4) 

Here, r^j is the distance between the z-th and the j-th monomer, 9{ and <pi are the 
bond and dihedral angles at the i-th monomer, respectively. The parameters a = 50 
and 7"o = 1 fix the strength of the harmonic force and the equilibrium distance 
between successive monomers. Both a and ro, as well as all the quantities in the 
following are expressed in dimensionless units, for a comparison with physical units 
see. 15 ) 

The value of a is chosen to ensure a value for V\ much larger than the other 
terms of potential in order to reproduce the stiffness of the protein backbone. The 
expression for the bond-angle potential term Vii^i) (|2-2p corresponds, up to the 
second order, to a harmonic term ~ kg(9{ — 9q) 2 /2, where 



cos(9 ) 

'sin 2 (0 o )' 4sin z (0 o ) 



A = -k 6 _ 2 K "°\ , B = - - ; /n - , V = Acas(p Q )+Boas(2e Q ) , (2-5) 



with kg = 20 and 9q = 57r/12 rad or 75°. 

The dihedral angle potential is characterized by three minima for p = (asso- 
ciated to a so-called trans state) and ip = ±2-7r/3 (corresponding to gauche states), 
this potential is mainly responsible for the formation of secondary structures. In par- 
ticular large values of the parameters Ci,D{ favor the formation of trans state and 
therefore of /3-sheets, while when gauche states prevail a-helices are formed. The 
parameters (Ci,Di) have been chosen in the following way: if two or more beads 
among the four defining <p are neutral (N) then Cj = and D{ = 0.2; in all the other 
cases Ci = Di = 1.2. The tapering function S(9i,9i + \) has been introduced in the 
expression of V3 in order to cure a well known problem in the dihedral potentials, for 
more details see. 7 )' 16 ) The quantity S{9i,9iJ r \) entering in the definition of V3 has a 
limited influence on the dynamics apart in proximity of some extreme cases. 

The last term V4, introduced to mimic effectively the interactions with the sol- 
vent, is a Lennard-Jones potential, which depends on the type of interacting residues 
as follows: if any of the two monomers is neutral the potential is repulsive cm,x = 
and its scale of energy is fixed by £n,x = 4; for interactions between hydrophobic 
residues cb,b = 1 and eb b = 4; for any polar-polar or polar-hydrophobic interaction 
cp 5 p = cp t B = —1 and £p t p = ep^p = 8/3. 

Accordingly, the Hamiltonian of the system reads 

h=k+v=y: vl ^ v Y vlt + e ^m+i) + 
1=1 t=i 

L-1 L-2 L-3 L 

i=2 i=2 i=l j=i+3 
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where all monomers are assumed to have the same unitary mass, and consequently 
the momenta can be defined as (p x ,i,Py,i,Pz,i) = (£i>yi,Zi)- 

In the present paper we consider the following sequence of 46 monomers: 

B g N 3 (PB) 4 N 3 B g N 3 (PB) 5 P 

This sequence that has been widely analyzed in the past for thermal fold- 
j n gi3)-i5),i7)-2i) ag we Y\ as f or mechanically induced unfolding and refolding. 7 )' 22 )' 23 ) 
The sequence studied exhibits a four stranded /3-barrel Native Configuration (NC), 
which is stabilized by the attractive hydrophobic interactions among the B residues 
(see configuration (a) in Fig. [3]). In particular the first and third Bg strands, forming 
the core of the NC, are parallel to each other and anti-parallel to the second and 
fourth strand, namely, (PB)± and (PB)$P. These latter strands are instead exposed 
towards the exterior due to the presence of polar residues. 

In the following we will report simulation results associated to two different 
kind of simulation protocols: equilibrium molecular dynamics (MD) canonical sim- 
ulations at temperature T performed by integrating the corresponding Langevin 
equation; steered out-of-equilibrium MD simulations intended to mimic the mechan- 
ical pulling at constant velocity of a protein attached to the cantilever of an atomic 
force microscope, or analogously when trapped in optical tweezers. In both cases 
the initial state of the system is taken equal to the native configuration (NC), that 
we assume to coincide with the minimal energy configuration. 



§3. Fluctuation relation and out-of-equilibrium unfolding 



Given a system with L particles, characterized by the Hamiltonian Ho(q), where 
q = {rj, pi} is a point in the system phase space, we are interested in evaluating the 
constrained free energy landscape 



Pfj(Q) = -In 



dq 5(Q - Q(q))e 



(3-1) 



where Q is some macroscopic observable, function of the microscopic coordinates q. 
If the system is driven out of equilibrium by an external potential U Z ^(Q), which 
depends explicitly on Q and on the external parameter z, whose temporal evolution 
is dictated by the protocol z = z(t), then the FEL (|3-ip can be obtained via the 
work fluctuation relation 5 -* 

(S(Q - Q(q))e-? w ) = e-4MQ)+^)(Q)]/Z . (3-2) 

In Eq. (|3-2p W is the work exerted on the system by the force associated with the 
potential U, and the symbol (•) refers to an average over all the possible stochastic 
trajectories spanning the system phase space, while the parameter z(t) changes over 
time. Zq is the partition function associated with the unperturbed Hamiltonian 
H (.q). 

Typically, the macroscopic observable characterizing the state of a biopolymer 
under mechanical stress is the end-to-end distance £ and this is the coordinate we will 
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consider in the following. Furthermore, we consider a quadratic potential U z n\ (£) = 
k/2(Q — z(t)) 2 , mimicking the effect of the force exerted on the molecule by the 
cantilever of an atomic force microscope, or by optical tweezers, where z(t) = v p x t, 
is the equilibrium position of the potential, moving with a constant velocity v p . 
We perform steered MD unfolding simulations, where one of the free ends of the 
molecules is kept fixed, while the other is pulled by the external force fc(£(t) — z{t)) 
associated with the potential U. Further details on the simulations are given in Ref. 7 ) 
We consider different values of v p , and for each of them we simulate a given number 
of unfolding trajectories. For each trajectory, we compute the work W done by the 
external force on the protein. Finally, in order to estimate fj(() from Eq. (|3-2j) we 
use the procedure introduced and discussed in Ref. 24 ) The results are plotted in 
Fig. [T]for different pulling velocities v p : we notice that as the velocity v p decreases, 
the curves collapse onto the same curve which corresponds to the best estimate of 
/j(C) given by the method described here. In Ref. 7 ) we have verified that the lower 
curve essentially coincides with an independent equilibrium estimate of the FEL 
obtained via the weighted histogram analysis method. 
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Fig. 1. Free-energy profiles fj as a function of the end-to-end distance £, for the model protein 
discussed in Sec. [2j as obtained by implementing Eq. (|3-2p . with T = 0.3, and for different 
pulling velocities: from top to bottom v p = 5 x fO -2 , f x f0 -2 , 5 x 10~ 3 , 5 x 10 -4 , 2 x f0~ 4 , 
2 x f0~ J and 5 x fO -6 . In (b) an enlargement of the curve for v p = 5 x fO -6 at small £ is 
reported. The number of different pulling trajectories considered to estimate the profiles ranges 
between 150 and 250 at the highest velocities to 28 at the lowest velocity v p — 5 x 10 -6 . The 
letters (a, b, c, d) indicate the value of /j(C) corresponding to typical configurations reported 
in Fig. [3] and the (blue) vertical solid lines the location of the Structural Transitions (STs), see 
discussion in Sec. [4] 



§4. Inherent structure formalism and structural transitions 

Inherent structures correspond to local minima of the potential energy, in par- 
ticular the phase space visited by the protein during its dynamical evolution can 
be decomposed into disjoint attraction basins, each corresponding to a distinct IS. 
Therefore, the canonical partition function can be expressed, within the IS formal- 
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Fig. 2. (Color online) (a) Free energy profile /js(C) as a function of the end-to-end distance, (b) 
average potential energy Vjs(C) vs C- The dashed vertical lines indicate the location of the three 
structural transitions. All the data refer to a data bank of ISs obtained via out-of-equilibrium 
steered MD simulations mimicking mechanical protein unfolding performed at T = 0.3 with 
v p — 5 x 1CF 4 . In the upper panel (a) it is reported for comparison also /j(C) (orange pluses) 
obtained via the fluctuation relation procedure with v p = 5 x 1CT 6 . 



ism, as a sum over the non-overlapping basins of attraction, each associated to a 
specific minimum (IS) a: 25 ^ 27 ^ 

Z IS (T) = J-^V^" f e-^ r Ur = Y^e-^ +R ^ (4-1) 

a ^ a a 

where N' is the number of degrees of freedom of the system, A is the thermal wave- 
length, r represents one of the possible conformations of the protein within the 
basin of attraction of a, V a is the potential energy associated to the minimum a, 
AV a (r) = V(r) — V a and R a (T) the vibrational free energy due to the fluctuations 
around the minimum. 

The free energy of the whole system at equilibrium is simply given by fis (T) = 
— Tln[Zis(T)]. However, in order to construct a free energy landscape as a function 
of a parameter characterizing the different IS, like e.g. the end-to-end distance £, 
it is necessary to define a partition function restricted to ISs with an end-to-end 
distance within a narrow interval [C; C + d£] 

Z IS (Q,T)=Y^^ l3[Va+Ra[T)] (4-2) 

a 

where the ^ indicates that the sum is not over the whole ensemble of ISs {a} 
but restricted. The free energy profile as a function of £ can be obtained by the 
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Fig. 3. (Color online) Typical configurations of the model protein along an unfolding trajectory 
driven by the mechanical force discussed in Section [3] with T = 0.3. The NC (a) has £o ~ 1-9; 
the other configurations are characterized by £ = 6.8 (b), £ = 16.8 (c), and £ = 27.1 (d). The 
beads of type N, B, and P are colored in green, red and yellow, respectively. 



usual relationship fis((,T) = —Tln[Zjs((,T)]; while the average potential energy, 
corresponding to ISs characterized by a certain £, can be estimated as follows: 

V IS ((,T )= ^ " Zis((T) . (4.3) 

In order to built a data bank containing the different ISs, we have performed 
mechanical unfolding simulations of the protein at different temperatures via steered 
Langevin MD integration schemes. The data bank contains 3, 000 — 50, 000 ISs 
depending on the examined temperature as detailed in Ref. 7 ^ It is worth to notice 
that the ISs have been collected by following out-of-equilibrium trajectories induced 
by mechanical manipulation on the NC at a velocity v p = 5 x 10 -4 . As shown in Fig. [2] 
(a) , the free energy profile reconstructed with this approach is almost coincident with 
the best estimate of fj(()- However, the fluctuation relation procedure, described 
in the previous Section, requires a pulling velocity which is two orders of magnitude 
smaller in order to obtain a reliable reconstruction. 

Moreover, referring to Fig.[T]and Fig. [2] (a), it is possible to identify the structural 
transitions (STs) induced by the pulling experiment. As shown in Fig. D](b), the free 
energy profile exhibits a clear minimum in correspondence of the end-to-end distance 
of the NC (namely, Co ~ 1.9). In more detail, up to £ ~ 5.6, the protein remains 
in native-like configurations characterized by a /3-barrel made up of 4 strands, while 
the escape from the native valley is signaled by the small dip at £ ~ 5.6 and it is 
indicated as ST1 in Fig. Q] (b) and Fig. [2] (a). 

For C > 6 the configurations are characterized by an almost intact core (made 
of 3 strands) plus a stretched tail corresponding to the pulled fourth strand (see 
configuration (b) in Fig. [3|). The second ST amounts to pull the strand (PB)^P 
out of the barrel leading to configurations similar to (c) reported in Fig. In the 
range 13 < ( < 18.5 the curve fis(0 appears as essentially flat, thus indicating 
that almost no work is needed to completely stretch the tail once detached from the 
barrel. The pulling of the third strand (that is part of the core of the NC) leads to 
a definitive destabilization of the /3-barrel. This transition is denoted as ST3 in Fig. 
[2] (a). The second plateau in fis(C) corresponds to protein structures made up of a 
single elongated strand (an example of this state is configuration (d) in Fig. [3j). 
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§5. Thermodynamical transition temperatures 

The main thermodynamic features of a protein can be summarized with refer- 
ence to three different transition temperatures: 7 )' 17 - 1 ' 20 )' 25 )' 28 ) the collapse temper- 
ature Tq discriminating between phases dominated by random-coil configurations 
rather than collapsed ones; the folding temperature Tf, below which the protein 
stays predominantly in the native valley; and the glassy temperature T g indicating 
the freezing of large conformational rearrangements. 27 ) Following the procedures re- 
ported in Ref., 28 ) we have determined these temperatures and obtained Tg = 0.65(1), 
Tf = 0.255(5), and T g = 0.12(2). These values are in good agreement with those 
reported in, 17 )' 20 ) where Tf and T g have been identified via different protocols. 

Then, we can try to put in correspondence the three unfolding stages previously 
discussed in Sec. 4 with thermodynamical aspects of the protein folding. In particu- 
lar, by considering the energy profile V/s(C) reported in Fig. [2](b), an energy barrier 
AVis and a typical transition temperature T t = (2AVis)/(3N), can be associated to 
each of the STs. The first transition ST1 corresponds to a barrier AVjs = 8(1) and 
therefore to Tt = 0.11(1), that, within error bars, essentially coincide with T g . For 
the ST2 transition to occur, the barrier to overcome is AVis = 16(1) and this is as- 
sociated to a temperature Tt = 0.23(2) (slightly smaller than Tf). The energetic cost 
to completely stretch the protein is 50(1) that corresponds to a transition tempera- 
ture T t = 0.72(1), that is not too far from the ^-temperature given above. At least 
for this specific sequence, our results indicate that the observed out-of-equilibrium 
STs induced by mechanical pulling can be put in direct relationship with equilibrium 
thermal transitions usually characterizing the folding/unfolding processes. 

§6. Dynamical transitions and structure fluctuations 

In the analysis of equilibrium properties of heteropolymers, an abrupt deviation 
of the structural mean square displacement from a linear temperature dependence is 
usually associated to a dynamical transition. 30 ^ Fluctuations of the protein structure 
at equilibrium have been studied experimentally via elastic incoherent neutron scat- 
tering as well as Mossbauer absorption spectroscopy. 30 )^ 32 ) These studies indicate 
the existence of different dynamical regimes separated by dynamical transitions, 30 ) 
in particular for hydrated proteins powders, a first non-linear enhancement of the 
mean square displacement with the temperature is observed around 150 K, and a 
second one around 240 K. The first dynamical transition is associated to torsional 
motions and observable also for dehydrated or solvent-vitrified system, while the sec- 
ond one is related to the onset of small-scale libration motions of side-chains induced 
by water at the protein surface. 

Fluctuations of the protein structure at a certain temperature can be character- 
ized in terms of the following indicator 27 ) 

1 L 

{(Au 2 )) = T Y. Au ^ where Au i = ( d lc M ) ~ &,cm) 2 , (6-1) 
1=1 
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associated to the fluctuations of the distance di,cM between the i-th residue and the 
center of mass of the protein, the symbols (•) refer to temporal averages, while ((•)) 
to an average over all the beads composing the heteropolymer. 

In the present case, we have estimated ((Au 2 )) for various temperatures T by 
performing equilibrium unfolding MD canonical simulations and by following the 
protein trajectory for a time t = 500, 000. In Ref. 27 ^ the authors have shown 
for off-lattice Go models, reproducing the B\ domain of protein G, the existence 
of a dynamical transition temperature Td ~ 0.4 x Tf denoting the onset of large 
scale fluctuations. As already mentioned in the previous Section, our model (as real 
proteins) is characterized by three different transition temperatures, at variance with 
the Go model examined in Ref., 27 ) where the folding and the collapse temperature 
coincide. It is therefore quite instructive to examine how many dynamical transitions 
are present in our model and their location in temperature. 

From Fig. U] (a) it is clear that ((Au 2 )) exhibits a linear behaviour until Tdi ~ 
0.2 ~ 0.78 x Tf where a sharp increase takes place. At T < Tdi the equilibrium 
dynamics is simply characterized by small harmonic oscillations of the beads around 
their equilibrium positions. Therefore by applying the theorem of equipartition of 
energy to the corresponding potential term (|2-ip we expect that 

((Au 2 )) = A£^i T = 7l T for T < Tdi , (6-2) 

as indeed verified (see Fig. IHa)). By approaching the folding temperature there is a 
strong nonlinear enhancement of ((Au 2 )) due to a configurational rearrangement of 
the protein structure, which can be associated to an activation process which leads 
the protein to cross the free energy barrier at ST2 (see Fig. [2ja)). 

Moreover just above Tf a second linear regime is observable. This is due to angu- 
lar oscillations around their equilibrium positions 6q, therefore by applying equipar- 
tition to the terms V\ and V2 in our model we derive the following dependence 

((Au 2 )) ~ 7l T + t^-^^T = (ti + 72 )T for T f <T < T D2 (6-3) 

and again this linear behaviour is in good agreement with the data up to Td2 ~ 0.5 ~ 
0.77 x Tq (see Fig. SJa)). In this temperature range, the protein is partially unfolded, 
it is no more in native-like configurations, and the degrees of freedom associated to 
bending fluctuations, involving three consecutive beads, are now activated. At the 
temperature Tdi we observe a second dynamical transition involving large configu- 
rational fluctuations. This dynamical transition can be considered as a precursor of 
the collapse transition, characterized by the complete unfolding of the protein and 
associated to the crossing of the free energy barrier at ST3 in Fig. [2ja). 

As a matter of fact just above Tq an almost linear regime is observable in a 
narrow temperature interval, namely 0.68 < T < 0.80, we believe that this further 
linear regime is due to fluctuations of the dihedral angles entering in the potential 
term V3. We found that, in this temperature range the linear increase is characterized 
by a slope (71 + 72 + 73) ~ 2.24, leading to an estimate for the new contribution 
73 ~ 2.07 which is of the order of 3/Di = 2.5 (by assuming that no neutral bead 
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Fig. 4. (Color online) Variance ((Au 2 )} versus temperature T, as estimated from equilibrium un- 
folding simulations of duration t — 500, 000. The (black) dots are the results of the simulations, 
while in (a) the dashed (orange) line corresponds to {{Au 2 }} = 71T with 71 = 0.02935 and the 
dash-dotted (magenta) line to {{Au 2 }} = (71 + 7 2 )T + d with 72 = 0.1434. In (b) the dotted 
(green) line refers to {{Au 2 }} — (71 + 72 + 7a)T + di with 73 = 2.067. 



is involved in the oscillating dihedral angles). These fluctuations are indeed more 
collective since they involve four consecutive beads. A third dynamical transition 
appears to take place around Tm ~ 0.82. This latter transition is probably related 
to fluctuations involving large part of the protein. 

Therefore each of the observed dynamical transitions is first characterized by 
small oscillations around some typical equilibrium configuration (this corresponds 
to the linear regime) followed by larger fluctuation induced by the breaking of hy- 
drophobic bonds and leading to a new equilibrium configuration of the protein (this 
phase is characterized by an abrupt increase in the protein fluctuations). Once the 
protein is rearranged one observes another linear regime due to the activation of a 
different set of degrees of freedom, which were previously hindered by the hydropho- 
bic interactions. At temperatures T < Toi the protein stays essentially in tightly 
packed native-like configuration and the only allowed oscillations are those of the 
beads around their equilibrium positions. At Tf < T < Td2 the protein visits a 
different sets of less packed equilibrium configurations of the free energy, which are 
still characterized by a native core essentially intact. In this regime the bending os- 
cillations of three consecutive beads become possible and are present together with 
harmonic oscillations of each bead. The transition at Tq leads essentially to config- 
urations almost completely stretched where fluctuations involving four consecutive 
beads (defining a dihedral angle) are now also activated. 

From the analysis of ((Au 2 )) we have no indication of the glassy transition 
occurring at T g , apart some fluctuation taking place just above T g as shown in 
Fig. (Ha). This is probably due to the fact that we have traced the dynamics of 
the protein for too short time windows. However, by increasing by a factor five 
the integration time, we do not observe substantial modifications in the behavior of 
((Au 2 )). 
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§7. Conclusions 

In the present paper, we have discussed how the FEL of a model protein driven 
out of equilibrium can be estimated by exploiting two different methods: namely, we 
applied a work fluctuation relation and the IS approach. The results obtained with 
the two methods compare well, although the IS approach provides a reliable estimate 
of the FEL already at larger pulling velocities compared to the first method. 

The FEL reveals three structural transitions along the unfolding pathways. By 
evaluating the potential energy landscape, we are able to assign a characteristic 
temperature to each of these structural transitions. Such temperatures compare 
well with the temperatures characterizing the thermal (un)folding of the molecule. 
Finally we analyze in detail the equilibrium structure fluctuations which mark the 
folding and collapse thermal transitions. Inspection of these fluctuations' variance 
allows us to identify the dynamical transitions which turn out to be precursors of 
the two corresponding structural transitions (namely, ST2 and ST3). 

In conclusion, our work provides strong evidence that, at least for the present 
protein model, the mechanical out-of-equilibrium unfolding pathways can be rec- 
onciled with the thermal folding and unfolding ones, provided that one performs 
a detailed analysis of the relevant quantities, namely the free and potential energy 
landscapes, and the thermal transition temperatures. 
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